Misleading signatures of quantum chaos 
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The main signature of chaos in a quantum system is provided by spectral statistical analysis of 
the nearest neighbor spacing distribution and the spectral rigidity given by As(L). It is shown 
that some standard unfolding procedures, like local unfolding and Gaussian broadening, lead to 
a spurious increase of the spectral rigidity that spoils the Aa(L) relationship with the regular or 
chaotic motion of the system. This effect can also be misinterpreted as Berry's saturation. 
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Quantum chaos has been an active research field since 
the link between energy level fluctuations and the chaotic 
or integrable properties of Hamiltonian systems was con- 
jectured 0J^], providing one of the fundamental signa- 
tures of quantum chaos (^^] in atoms, molecules, nuclei, 
quantum dots, etc. The secular or smooth behavior of 
of the level density is a characteristic of each quantum 
system, while the fluctuations relative to this smooth be- 
havior are related to the regular or chaotic character of 
the motion in all quantum systems. To achieve the sepa- 
ration of the smooth and fluctuating parts, the energy 
spectrum is scaled to a sequence with the same local 
mean spacing along the whole spectrum. This scaling 
is called unfolding |5j . Although this can be a non-trivial 
task H , the description of the unfolding details of calcu- 
lations is usually neglected in the literature. 

In this Letter we show that, contrary to common as- 
sumptions, the statistics that measure long range level 
correlations are strongly dependent on the unfolding pro- 
cedure utilized, and some standard unfolding methods 
give very misleading results in regard to the chaoticity 
of quantum systems. Long range level correlations are 
usually measured by means of the Dyson and Mehta A3 
statistic ||. On the other hand, short range correla- 
tions, caracterized by the the nearest-neighbor spacing 
distribution P(s), are not very sensitive to the unfolding 
method. 

Let us consider a rectangular quantum billiard with a 
size ratio a/b = n. This is a well known example of a 
regular system. In general, for regular systems level fluc- 
tuations behave like in a sequence of uncorrelated energy 
levels, and the A^(L) statistic increases linearly with L. 
However, it was shown by Berry that the existence of 
periodic orbits in the phase space of the analogous clas- 
sical system leads to a saturation of A^(L) for L larger 
than a certain value L s , related to the period of the short- 
est periodic orbit. Fig. |l| shows the A3 behavior for 
a sequence of 8000 high energy levels of the mentioned 
quantum billiard, calculated with two different unfolding 
procedures. The mean level density for this system is 
given by the Weyl law Q . Using this density to perform 
the unfolding, A 3 follows the straight line of level spac- 



ings with Poisson distribution, characteristic of regular 
systems. In this example, Berry's saturation takes place 
at L s ~ 750, that is outside the figure. Let us suppose 
now that the law giving the mean level density of the sys- 
tem were unknown. Then, a standard method to obtain 
the local mean level density at energy E is to calculate 
the average density of a few levels around this energy. 
Using this method one obtains a very different behavior, 
the spectral rigidity increases strongly at L ~ 20, and af- 
terwards A3 is close to the Gaussian orthogonal ensemble 
(GOE) line characteristic of chaotic systems. The latter 
behavior is not at all related to the Berry saturation. 
It is a spurious effect due to inappropriate unfolding of 
the level spectrum and it implies that strong long-range 
correlations have been improperly introduced by the pro- 
cedure. 
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FIG. 1. Comparison of the A3 statistic for a rectangular 
quantum billiard using two unfolding procedures. Open tri- 
angles correspond to the smooth unfolding using the Weyl 
law, and filled triangles to the local unfolding method. The 
dashed line is the Poisson limit and the dotted line is the GOE 
limit. 

This first example illustrates the problem that can 
arise with some reasonable unfolding methods currently 
used in quantum chaos calculations p|-|l3|]. In order to 
understand its origin we shall analyze different unfold- 
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ing procedures. The principal difficulty in the unfolding 
is the correct characterization of the mean level density 
p(E). Having this function, the unfolded adimensional 
variables Si, 



£l = N{Ei), N(E) 



dEp(E), 



(1) 



have mean level density p(e) = 1. The unfolded spac- 
ing sequence is then {.Si — £i+\ — and the nearest- 
neighbor spacing distribution P(s) is well suited to study 
the short-range spectral correlations Q . 

The A3 statistic is used to investigate the long range 
correlations. It is defined for the interval [a, a + L] in the 
cumulative level density as 



A 3 (a,L) = - 



mm 

A.B 



[N(e) - As - B} 2 de. (2) 



The function A^^L), averaged over intervals, measures 
the deviations of the quasi-uniform spectrum from a true 
equidistant spectrum. 

For some systems a natural unfolding procedure ex- 
ists, because p(E) is known from an appropriate statis- 
tical theory or by a well checked empirical ansatz. For 
example, p(E) is a semicircle for large GOE matrices || , 
it often has Gaussian form for large nuclear shell-model 
matrices ||, and follows the Weyl law in quantum bil- 
liards ||. However, in many systems where there is no 
natural choice for p{E), it is usually estimated from a 
set of neighboring levels. The simplest method, called 
local unfolding, has been widely used [p[-^2[. The mean 
level density is assumed to be approximately linear in a 
window of v levels on each side of Ei , and is given by 



Pl{E % ) = 



2v 



Ei-\ 



(3) 



where L stands for local unfolding. More sophisticated is 
the Gaussian broadening method ]8|,|T3|]. The level den- 
sity p{E) = 3(E ~ Ei) i s substituted by an average 
level density 



Pg{E) = 



1 



5>xpi - 



( E-Etf 
2a 2 



(4) 



where G stands for Gaussian broadening. The sum 
runs over all the energy levels, but only those satisfy- 
ing \E — Ei\ < a do significantly contribute to p@(E). 
Although these two methods are different, both depend 
on a parameter v or a that measures, in a real or effective 
way, how many neighboring levels are used to calculate 
the local mean density. 
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FIG. 2. A3 for the complete J = 10 level sequence of a 
shell- model calculation for 52 Ca in the pf shell. For smooth 
unfolding made by an Edgeworth expansion in the cumulants, 
the result (open triangles) lies between Poisson (dashed line) 
and GOE (dotted line) limits. Filled triangles correspond to 
local unfolding with a 2v = 10 window. 

Let us now consider the atomic nucleus as example of 
a quantum system more complex than the quantum bil- 
liard. In most nuclei, level fluctuations are in agreement 
with GOE predictions at all energies, showing that the 
motion is chaotic. However, it has recently been observed 
that single closed nuclei are less chaotic than expected 
0,[L4|. One of the most regular nuclei at law energy 
is> ^7a. Analysis of the shell-model level spectrum |l4| ] 
shows that the nearest neighbor P(s) distribution is close 
to the Poisson limit (the Brody parameter is u> = 0.25) 
for levels up to 5 MeV above the yrast line. As the exci- 
tation energy is increased, P(s) approaches the spacing 
distribution of a chaotic system. However, other statis- 
tics indicate that the dynamics still is not fully chaotic. 
Often, the mean level density inside a valence space is 
very well reproduced by an Edgeworth expansion around 
a Gaussian form JsJ]. If we use it to perform a smooth 
unfolding of the 2755 J — 10 levels of 52 Ca, the A 3 (L) 
statistic is close to GOE limit for very small L values, but 
it increases linearly instead of logarithmically for larger L 
values, as can be seen in Fig. ^. Except for very small L 
values, the spectral rigidity is intermediate between that 
of GOE and Poisson limits, giving a clear signature of 
non chaotic motion. This result is in agreement with the 
behavior of the wave function localization lengths [O] . 

When the Edgeworth expansion fails, as it happens 
sometimes p| , the local unfolding or the Gaussian 
broadening are the available unfolding methods. Fig. |^ 
also shows the results of local unfolding for v = 5. The 
calculated A3 follows the line obtained with the smooth 
unfolding up to L ~ 2v, but then accumulated unfold- 
ing errors increase the spectral rigidity and lead to a A3 
saturation for larger L. This is the same behavior that 
was observed in the quantum billiard system described 
above. Moreover, as the A 3 values are rather close to 
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the GOE limit when v = 5 is used, the conclusion in this 
case would be that 52 Ca is a chaotic system. This ex- 
ample is then very enlightening. First, it illustrates that, 
contrary to common practice, the A^(L) statistic should 
be calculated up to high L values, because otherwise one 
can miss relevant information on the system dynamics. 
Second, it shows the problems that can arise when the 
mean level density is not known and one has to rely on 
local unfolding. 

To avoid any uncertainties on the real mean level den- 
sity p(E) and level fluctuations, we can study GOE and 
Poisson level spectra, the paradigmatic cases of chaotic 
and regular systems, respectively. We consider a GOE 
matrix with dimension N = 10000, and compare the 
spectral fluctuations obtained by three different meth- 
ods: Smooth unfolding made with the semicircle law, 
local unfolding, and Gaussian broadening unfolding. All 
these methods yield almost indistinguishable results for 
the P{s) distribution, that is in perfect agreement with 
the Wigner surmise. The behavior of the short range 
correlations is not affected by the method of unfolding. 
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FIG. 3. Comparison of the spectral rigidity for an 
N — 10000 GOE level spectrum, calculated with different 
unfolding methods. Open triangles correspond to the smooth 
unfolding. Filled circles and squares to the local unfolding 
with window size 2v = 42 and 2v = 18, respectively. Open 
circles and squares to the Gaussian broadening unfolding for 
cr = 2 and a = 1 MeV, respectively. 

However, Fig. ^ shows a completely different scenario 
for A3. For the smooth unfolding, the spectral rigidity 
behaves as predicted by GOE, up to very large L values. 
The local unfolding was performed using two different 
windows, with v = 9 and v = 21. The calculated A3 
coincides now with GOE predictions only up to L ~ 2v, 
then it leaves the GOE trend because the spectral rigid- 
ity increases and finally A 3 saturates to a constant value. 
The Gaussian broadening unfolding was performed for 
(7 = 1 MeV and a — 2 Mev. In the central part of the 
spectrum these values correspond to windows contain- 
ing about 10 and 20 states, respectively . Therefore, the 



effective number of states that affect the average level 
density is about the same as in the local unfolding case. 
Again, we see the same A3 behavior for L values greater 
than the window used in the unfolding. 

Fig. [| shows the spectral rigidity for 10000 levels 
generated with Poisson statistics and a uniform density 
p(E) = 1. The smooth unfolding gives A3 values close 
to Poisson predictions, but local unfolding with v = 2, 
9 and 21, leads again to the same behavior observed in 
previous cases for L > 2v. In fact, for the small win- 
dow with four spacings, the A3 of the Poisson spectrum 
closely follows GOE predictions! 
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FIG. 4. Comparison of the spectral rigidity for a Poisson se- 
quence of 10000 levels with uniform density, using several un- 
folding procedures. Open triangles correspond to the smooth 
unfolding. Filled circles, squares and triangles, to the local 
unfolding with window size 2v — 42, 2v = 18 and 2v — 4, 
respectively. 

Looking for deeper insight into the spurious A3 sat- 
uration, we consider the sequence of nearest level spac- 
ings as a physical signal, and apply Fourier analysis tech- 
niques to its study. We have chosen a system with Pois- 
son statistics and uniform level distribution to illustrate 
the idea, because the smooth density is constant. There- 
fore, the fundamental assumption of the local unfolding 
method, namely that the mean density is approximately 
linear within a window, is exactly fulfilled. From the real 
nearest-neighbor level spacing sequence S, we obtain: (a) 
the average spacing sequence Dl calculated with the lo- 
cal constant density of Eq. (g) and v — 21, (b) the 
sequence of smoothly unfolded spacings s, and (c) the 
sequence of locally unfolded spacings sl- Since for this 
spectrum p(E) — 1, we have D — 1 and s — S. 

Fig. displays the power spectrum of these sequences 
for frequencies up to k = 0.6. For Dl, it has a maximum 
near k = and decreases smoothly becoming essentially 
zero at some threshold frequency ko — ir/v. However, 
this behavior is a spurious effect, because the real mean 
spacing D is constant and then its power spectrum is 
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zero for all the frequencies k ^ 0. Therefore the local 
unfolding procedure introduces spurious low frequency 
components into the D L signal. Comparing the power 
spectra of Dl and s, it is seen that they are very similar 
at low frequencies, except for the damping of the former. 
The power spectra of s and sl are also very similar, ex- 
cept that the low frequency components are missing in 
the latter. These results clarify the deficiencies of local 
unfolding. It becomes apparent that the procedure is fil- 
tering out low frequency fluctuations from the spectrum 
s, and improperly including them in Dl- Moreover, by 
reducing or eliminating fluctuations of frequency smaller 
than ho, the procedure is introducing long range corre- 
lations with wave lengths greater than 2v. As this fluc- 
tuation reduction is progressive, the spurious long range 
correlations become stronger as L increases beyond the 
window size 2v. It is precisely this phenomenon what has 
previously been detected by the A3 statistic, that strong 
long range correlations leading to a saturation of the A3 
are observed for L > 2v. 



small windows the behavior of A3 may be close to the 
GOE limit. Furthermore, the spurious saturation of A3 
can easily be misinterpreted as Berry's saturation. 

For systems intermediate between regular and chaotic, 
the traditional spectral statistics P(s) and As(L) for 
small L values may be close to the GOE limit, and 
strong deviations of the spectral rigidity from GOE pre- 
dictions only appear for larger L values. Thus if the 
local mean level density is not known from a statistical 
theory or a good empirical ansatz, the analysis of energy 
level fluctuations will not lead to correct conclusions on 
the system dynamics. In this case, it becomes neces- 
sary to go beyond level statistics and study properties of 
the wave functions, such as localization length, transition 
strengths and transition strength sums |l5|] . 
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FIG. 5. Comparison of the power spectrum for the se- 
quences Dl, s and sl, for a sequence of 10000 uncorrelated 
levels with uniform density. The locally unfolded spacings sl 
are calculated with a window of size 2v — 42 that corresponds 
to fe = 0.15. 

In summary, we have shown that the correct behav- 
ior of A3 is strongly modified by some commonly used 
unfolding procedures when the exact shape of the mean 
level density is not known. Methods like local unfolding 
or Gaussian broadening introduce spurious long range 
correlations in the unfolded level spectrum, increasing 
the spectral rigidity and leading to a saturation of A3(i). 
In these methods the local average level spacing at energy 
E is calculated from the levels inside an energy window 
around E. The spurious behavior of the A3 statistic is 
observed for L larger than the window size. In general 
it gives misleading signatures of quantum chaos, and for 
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